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We introduce the quantum mechanical formalism for treating surface plasmon polariton scattering at an in- 
terface. Our developed theory - which is fundamentally different from the analogous photonic scenario - is used 
to investigate the possibility of plasmonic beamsplitters at the quantum level. Remarkably, we find that a wide- 
range of splitting ratios can be reached. As an application, we characterize a 50:50 plasmonic beamsplitter and 
investigate first-order quantum interference of surface plasmon polaritons. The results of this theoretical study 
show that surface plasmon beamsplitters are able to reliably and efficiently operate in the quantum domain. 
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I. INTRODUCTION 

Nanophotonic systems based on surface plasmon polari- 
tons (SPPs) |Q1 0] are currently raising considerable interest 
from the quantum optics and quantum information commu- 
nities O- ITori . Due to their ti ght- field confinement lUlll and 
electro-optical behavior yj, |2[ [l2|], SPPs constitute compact 
and versatile candidates for quantum information processing 
(QIP) with light at the nanoscale. One of the most impor- 
tant ingredients for plasmonic -based QIP is the ability for 
SPPs to interact coherently with each other. Recent theoreti- 
cal Il3l 1 1 4f1 and experimental work 1 1 51 - 4 1 811 has hinted at the 
possibility of achieving coherent interactions between SPPs 
via scatte ring ty pe processes. Here, metal-dielectric inter- 
faces Il3lll4ll9l] . as well as junctions and splitters in waveg- 
uides structures il7l[l8ll have been considered. However, this 
work has so far been restricted to a purely classical domain 
and little is known about surface plasmon interactions at the 
quantum level. The development of a flexible theory that can 
be applied to a variety of different types of waveguide geome- 
tries should greatly aid the design of plasmonic components 
exploiting quantum mechanical effects and circuitry for QIP 
applications. 

In this work we introduce and develop the quantum me- 
chanical formalism for treating SPP interactions via scatter- 
ing at an interface. Our theory is then used to investigate the 
possibility of constructing plasmonic beamsplitters that are 
able to operate faithfully and efficiently at the quantum level. 
Surprisingly, we find that a wide-range of splitting ratios can 
be reached without the need for embedding complex optical 
material, such as anisotropic metamaterials, as recently sug- 
gested [14]. We also find that power loss due to unavoid- 
able so-called 'parasitic scattering' of SPPs into photon ra- 
diation [13] can be suppressed to 5% or even less in some 
cases. Moreover, the beamsplitter geometries we investigate 
are directly accessible to experiments; due to the basic proper- 
ties of the materials involved, complicated on-chip fabrication 
techniques are not required. As an application of our theory 
- and as an example of the necessity for a quantum theory of 
scattering - we optimize a 50:50 plasmonic beamsplitter and 
use it to investigate first-order quantum interference effects of 
SPPs. Our study shows that surface plasmon beamsplitters 



can reliably operate at the quantum level, providing important 
insights and helping to open up new directions of research into 
the design of efficient and practical components for on-chip 
plasmonic-based QIP. 

The paper is structured as follows: In Section II we intro- 
duce the interface scenario considered for SPP interactions via 
scattering. Here, a quantized description of all the fields in- 
volved is presented and a brief discussion of relevant material 
properties is included. In Sections III and IV we formalize our 
theory for the scattering process, where we provide a normal- 
ization procedure for the quantized fields and introduce field- 
matching relations. This enables the formation of a quantum 
transfer matrix in Section V linking all the relevant excita- 
tions together. In Section VI we use our theory to investigate 
the possibility of efficient quantum plasmonic beamsplitters. 
Here, we discuss important issues such as spatio-temporal in- 
distinguishability, reciprocity and loss effects. In Section VII 
we use the results from the previous sections to investigate 
first-order quantum interference of SPPs. Finally, Section VIII 
summarizes our main results. 



II. INTERFACE CONFIGURATION 

The interface considered is shown in Fig.[T](a). Here two 
regions i and j, with different materials are combined at x = 
(inset shows the cross-section). Region i (J) consists of metal 
with permittivity e m j (e m j) for z < and a dielectric media 
with permittivity e^, (edj) for z > 0. Before we introduce the 
quantum formalism, a short description of the system dynam- 
ics is given. The purpose of this is to provide an informative 
glimpse of the detailed theory that is to follow. 

In Fig. Q](a) an SPP excitation, denoted by the operator a?, 
is shown in region i moving forward from the left at an angle 
6> 1; . This SPP scatters at the interface (x = 0) into a backward 
moving SPP excitation, denoted by the operator a b , in region 
i at an angle 6>, r , in addition to a forward moving SPP excita- 
tion, denoted by the operator V , in region j at an angle 0,- f , 
as well as into near- and far-field photon radiation excitations 
moving backwards and forwards in regions i and j respec- 
tively, denoted by the operators A b and W (all angles are in 
the x-y plane and shown in more detail in Fig.|2](b)). The ra- 
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FIG. 1: (Color online) Scattering of SPPs at an interface, (a): An SPP excitation a! in region i is scattered into an SPP excitation V in region j 
and SPP excitation a b in region i. The superscript / (b) denotes a forward (backward) propagating excitation with respect to the x axis. Near- 
and far- field radiation excitations are also excited in the process (represented by the central white haze and red jagged arrows). The symbols 
A b k and are used for the radiation in regions i and j respectively. Similar scattering occurs for the SPP excitation b b in region /, Further 
details on the angles are given in Fig. 2 (b) and related discussion in Section V. (b): Dispersion relations for the SPP's and radiation in regions 
i and j. Shaded regions correspond to a continuum spanned by the wavenumber q\ and qj respectively. Examples of lines corresponding to set 
qi and qj are shown. 



diation excitations, have a range of in-plane and out-of-plane 
wavevector components (with respect to the x-y plane) and 
act as a loss mechanism for SPP interactions at the interface. 
Later, we will show how this loss can be suppressed to less 
than 5% by carefully modifying the properties of the metal 
and media in regions i and j, as well as the incidence angle of 
the incoming SPP. A similar outline to the above can also be 
given for an SPP excitation moving from region j to i. 

We now formalize the above scenario. In order to describe 
the scattering process quantum mechanically, we start by in- 
troducing the quantized vector potential field operators for the 
various excitations supported by the geometry of Fig. [T] (a) 
and involved in the interaction. We consider the vector po- 
tential as it provides a convenient starting point from which 
to derive both the quantized electic and magnetic field opera- 
tors l20ll . These operators will play a central role in our theory. 



A. Quantized radiation field 

In region i both transverse magnetic (TM) and transverse 
electric (TE) surface radiation fields are supported I2TI l22tl . 
This is in contrast to SPPs which are restricted to TM fields, 
as we will discuss in the next section. The quantized vector 
potential for the radiation fields is derived from the field equa- 
tion for the geometry in the Coulomb gauge with classical- 
quantum correspondence relations applied 12111 . It is given by 



A„(r, t) 



'twfqj 



£Q7TCO 



(1) 



[0 r (k ; , q h r)A r (k,-, <?,>-""<' + #(k,, q h r)A r T (k ; , qi)e m>, l 

The creation (annihilation) operators A,t(k, qi) (A r (k,-, qi)), 
with r € {tm, TE}, satisfy the bosonic commutation relations 



[Arik^qilA^q'i)] = 6(ki-tyd(qi-q'Jd rr > from which the 
Heisenberg uncertainty relations are derived [20]. The wave- 
functions are given by 



0™(k ! -, 9i ,r) = e * r r [ M 



(ifci + -z)e ViZ §{-z) + 

Vi 
k; 



(2) 



[/k,(cos qiz - rjj sin q t z) + — z(sin q ( z + q, cos <7,z)]#(z) 
qi 



and 



n ,(ki,9i,r) = «*- r of »(* x ki)[e ViZ &(-z) + 

Vi 

(cos qiz + — sm qiz)v(z)], 



(3) 



with #(z) representing the heaviside function, k; = k x jx + kyjf 
is a wavevector with wavenumber k, = {k 2 xj + fc^) , qi = 
(utfedj/c 2 — k 2 ) 1 / 2 characterizes the variation of the wave- 
functions in the region z > and v, = (k 2 — a> 2 e,„j/c 2 ) 1 ^ 2 
characterizes the z < variation. In addition, e^, is the di- 
electric function of the non-metallic media in region i and 
e m j = 1 — (iXpJu 2 is the dielectric function of the metal, 
with a> Pi as the plasma frequency in region i and the pa- 
rameters T]i = q i e m ,i/v i £d,i, 7™ = ViKe^ .v 2 + e 2 m -q 2 ) 112 and 

jf = [e d M,i ~ £m,dr m JHHH. Here, the Drude model H 
is used in e mji for the purpose of illustrating the results of our 
investigation and we have extended the work of Ref. lf2lll to 
allow for an arbitrary dielectric function Q , that is real and 
positive. The parameter qi is chosen to be the free parameter 
characterizing the complete z variation of the wavefunctions. 
This can be done by substituting k; in terms of q-, in the def- 
inition of Vi giving vi = {u) 2 (e d j - e m ,i)/c 2 - q 2 ) 1 ' 2 - Thus 
the definition of g, provides a dispersion relation (DR) which 
links oji to k, regardless of the polarization (TM or TE). This 
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DR is set by the value of q,, once edj is specified. In Fig.Q](b) 
we show example DRs of the radiation in separate regions i 
and j for edj = 3 and edj - 1 . For q — 0, the radiation has the 
usual bulk photon DR expected from a medium with dielectric 
function ej, given by k = y/edoi/c. However, as q increases, 
the penetration of the field into the metal, characterized by the 
e vz term in Eqs. (f2]l_and ([3]), increases and the DR is modi- 
fied substantially 02111 . The shaded regions denoted by rad,- 
and rady correspond to a continuum of surface photon DRs 
spanned by q t and qj respectively. 

We limit our work to frequencies below oj p j. Here the inte- 
gral over qi in Eq. ([TJ explicitly covers the range to a max- 
imum cutoff value q cut = U)i(edj) 1 ' 2 /c. For q t > q cut , ki be- 
comes imaginary and for qi > cjjfaj — e m j) l ' 2 /c > q cut , ki 
and v, are imaginary. Both these ranges require a different an- 
alytical form for A,,(r, t). For frequencies above u> p j we have 
9cut > t (edj - <w) 1/2 and the regime ^(e d j - e„,,,) 1/2 < qi < 
9cut (where ki is real but v* is imaginary) corresponds to the 
transparency region of the metal, as shown in Fig.[T](b). The 
quantized vector potential for the radiation field in region j is 
given by Eqs. (Q]), (f2| and (f3]l with the relabeling i — » j and 
the creation (annihilation) operators A\ — > b\ (A r — > B r ). 



B. Quantized SPP field 

The quantized vector potential for the SPP field in region i 
is derived in a similar way to the radiation field IT7l[8L l2l[l24ll . 
However unlike the radiation field, due to the materials and 
geometry considered here, SPP boundary conditions lead only 
to TM fields being supported ujj, Ll3|, Ll4f| . The quantized vector 
potential is given by [21] 



A p ,(r, t) = 




with creation (annihilation) operators a^(k,) (S(k,)) satisfying 
the bosonic commutation relations [a(k,), a^(k,)] = 5(k, — k-) 
and wavefunctions 

</v(k,-,r) = [(ik f - + -z)e v ' z m-z) + (5) 

V; 

(i% - — z)e- V0 ^(z)]e*'- r . 

Vo,i 

Here, k, = k x jx + kyjf is again the wavevector with wavenum- 
ber ki = (k\ . + k%j) m , v ,i = (*• - ajre d j/c 2 ) 112 character- 
izes the variation of the wavefunctions in the region z > 
and V,- = (k 2 — afiemj/c 2 ) 1 ' 2 characterizes the z < varia- 
tion. The parameter p t = cti[2coief ni ed,i(-(e„,j + edj)) 1 ' 2 ]' 1 , 
with ti = (e 2 n i + €d.i)(e 2 ■ - e%J). Again, here we have ex- 
tended the work in Ref. i2lll in order to allow for an arbitrary 
real and positive q , . A DR that is independent from vo,; and 
V, can be derived from boundary conditions and is given by 
ki = (^i/cXedje^i/iedj + e m ,,)) 1/2 JHHMS. Once k t is 
set, with edj and e m j both specified, w, is set and so are v; and 



vo,, by definition. This is in direct contrast to the radiation ex- 
citations, where for a set wavenumber ki there are a range of 
values of which satisfy the DR, dependent on the free pa- 
rameter qi. Here, as v,- and vo,, are both set by ki, for SPPs there 
is no longer a free parameter for the z variation. Thus neither 
v, nor vo,, appears explicitly as an argument of the wavefunc- 
tions. In Fig.Q](b) we show example DRs of SPPs in regions 
i and j for edj = 3 and e^j = 1, with e m chosen to be that of 
silver, co pJ = co pJ = 1.402 x 10 16 rads~' JH. The quan- 
tized vector potential for the SPP field in region j is given by 
Eqs. © and (0 with the relabeling i — > j and the creation 
(annihilation) operators a* — > b* (a — > b). 

By inspection, the wavefunctions for the radiation and SPP 
excitations for a given region i are found to be orthogonal J2lll 
and form a complete set of eigenfunctions for the physical 
space and parameter regime considered B26I1 . Both these prop- 
erties will be essential for obtaining a consistent map be- 
tween the fields in two different regions. As can be seen 
from Fig.Q](b), in a given region the radiation and SPP DRs 
never cross; a manifestation of the orthogonality of the under- 
lying wavefunctions. Thus mode-matching conditions cannot 
be met. As a result, it is not possible for SPP and radiation 
excitations to directly couple to each other in the same re- 
gion. However, the situation changes with the introduction of 
an interface where two regions with different physical prop- 
erties are joined together. Depending on the material prop- 
erties, field-matching and energy conservation, coupling be- 
tween SPPs and radiation within the same and across two dif- 
ferent regions can occur. Before we treat such a situation in 
detail, we first ensure that all the excitations involved (as de- 
fined by Eqs. (Q]) and ©) are normalized correctly. 

III. NORMALIZATION 

A standard approach in classical coupled mode theory is to 
use the electromagnetic reciprocity and Poynting's theorems 
to ensure that the various fields are normalized correctly with 
respect to energy transfer in the direction normal to an inter- 
face i22ll . In the quantum case we check the normalization 
by starting from the quantized version of the Poynting vec- 
tor i20ll . which remains valid for the parameter regime consid- 
ered here, corresponding to a weakly dispersive media lflll27l- 
[30h - The quantized Poynting vector is given by 

S^(r, t) = &-(r, t) x H+(r, t) - ft-(r, t) x fe+(r, t). (6) 

Here, /u £ {r,p} corresponds to radiation (r G {tm, TE}) 
or SPP (p) excitations. The quantized fields E p i(r, t) and 
H^r, t) can be obtained in the usual way from the vector po- 
tential A^j(r, f) defined in Eqs. ([1) and <j4j with the relations 
E^,(r,f) = -d,k pi (r,t) andft^r.f) = ^'(V x A pi (r,t)). 
The same can be carried out for the quantized fields in region 
j. The ± superscript corresponds to the positive and negative 
frequency parts of a given field operator, the explicit form of 
which are provided in Appendix A. Upon substitution of the 
relevant fields into Eq. <JSJ and taking the x component (di- 
rection normal to the interface) integrated over time and y-z 
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cross-section at the point x — 0, one finds 
JJJ S^,-(r,0 xdtdydz 

= wjJ dMk ''Jo Jo «jJJ d ^ 
(E:(k i , ?i ,r,0xH;.(k^;,r,() 

-ft-.(k f , r, f) x E;,.(k;, r, t)) ■ x, (7) 

where for the SPP excitations (jx = p) the parameter q, and its 
integration are removed from the field definitions Dili . The 
limits on all integrals are — oo — > +oo unless stated otherwise. 

The reciprocity theorem of electrodynamics imposes an or- 
thogonality relation for the wavefunctions of the electromag- 
netic excitations with the same time dependent harmonic evo- 
lution l22ll . If we consider that the radiation and SPP excita- 
tions are not damped in their direction of propagation, i.e. k x j 
and k y j are both real (we will return to this point in more detail 
in Section VI), the orthogonality condition is given by 

fff (k,., q h r, t) x H v + ,.(k;, q[, r, t) ■ x dt&ydz 

= |A/;,(k ; , ft)| 2 c5(k - k^ft - q 1 ^, (8) 

where |A^,,(k/, qd\ 2 = 2n 2 ha>i is found using the relations pro- 
vided in Appendix A [32]. This expression allows us to recast 
Eq. (0 for forward propagating excitations (k x /. +ve) as 



III 



S„(r, t) ■ x dfdydz 



JJo ' 



dk,dft tkOi n r (k, ft), (9) 



for the radiation and 

Jjj S pi (r, t) • x dfdydz = J dk, Tkj; n^k,), (10) 

for the SPPs. Here h p = cfia and h r = A}A r are the par- 
ticle number operators. Physically, the right hand sides of 
Eqs. (|9]l and ( fTOb represent the sum of the excitation energies 
of the forward propagating field modes above their zero-point 
value [29]. Thus they represent the total radiation and SPP 
energy that flows through the y-z plane in the +x direction at 
x — 0. Similar equations hold for the backward propagating 
excitations. With the fields normalized correctly with respect 
to energy flow across the interface, we can now proceed to de- 
velop the formalism for the coupling between the excitations. 



and requires a great deal of care to be taken in the derivation 
due to inter-polarization coupling. 

For simplicity, the wavefunctions of all excitations are 
taken as plane waves, extending infinitely in the direction per- 
pendicular to that of propagation. A beam-width can be im- 
posed in a straightforward manner in order to make them finite 
and closer in description to an experiment, however this does 
not result in any change to the underlying theory B7ll8l l20ll33n . 

The total normalized quantized electric field on side ; con- 
sisting of SPPs and radiation can be written as 



E,(r,0 



[E;/'(r,k,f) + E/(r,k,f)]dk 



kyi+ve 
k r! +ve 



+ | [E;f(r,k,f) + E-. fe (r,k,f)]dk 



J [E;,°(r,k,f) + E- (r,k,f)]dk- 



ky • — ve 



nffv 



(r, k, ft, f) + E n /(r, k, ft, f)] dkdft 



JOT 



+ I I [E+V, k, ft, t) + E" D (r, k, ft, f)] dk,dft 



II 



+ [E+°(r, k, ft, t) + E-°(r, k, ft, t)} dk,dft 



(11) 



kyj — ve 



where the k, wavevector integral has been split into three parts 
for the SPPs and the radiation (TM and TE). The range of the 
components of the wavevectors for each of these parts is spec- 
ified below the integral sign. Here and elsewhere the super- 
scripts / and b on the operators denote 'forward' and 'back- 
ward' propagating excitations, whereas o denotes excitations 
with all 'other' propagation directions. As a special case, the 
superscript b denotes a reversing of the k x component in an 
operator and its associated wavefunction. This allows the inte- 
grals for the forward and backward excitations to be combined 
later. No change to operators or wavefunctions for the / and o 
excitations should be made. For the quantized magnetic field 
we have 



IV. QUANTIZED FIELD MATCHING 

We start by using the E and fields to derive the cou- 
pling between SPP and radiation excitations on either side of 
the interface shown in Fig. Q] (a). In doing this, we extend 
a recent classical study on coupled hybrid modes lfT3ll to the 
quantum domain. At the same time, the theory we develop 
here is more general than Ref. [ 13] in the sense that it includes 
SPP scattering at an angle of incidence to the interface. This 
is not a straightforward extension, even in the classical case 



[ft^Cr.ki.O + ft/Cr.ki.Oldki 



H,(r,?)= J 

I 

>v' +ve 

c j +ve 

+ J [H;°(r,k,f) + V (r ' k '' f)]dki 



f f [H;f(r,k,f) + H- 6(r ' k '' f)]dki 

it r ; +ve 



5 



+ J] [ jf [ftV (r, k„ ft, f) + H/(r, k„ ft, f)] dMft 
+ Jj [H^(r, k„ ft, f) + ft-»(r, k„ ft, f)] dkjd^j 

ft ,'+ve 

+ Jj Dft^(r,k,, 9 i,0+ftrt O (r,ki,?i,0] dk,d<? ; 



(12) 



ky j — ve 
* ';±ve 



Similar expressions can be written for the quantized fields on 
side j. We now match the transverse part of the fields on either 
side at point x = as follows 



<«S(r,f) = ^(r,f), 
Mr,t) = J^(r,/), 



(13) 
(14) 



where $(r, f) = (E,(r, t)-f)f +iE,(r, f) • z) z and ^(r, f) = 
(H;(r, t) • y)y + (H,(r, f) • z)z [34]. By noting the following 
relations for the transverse component of the wavefunctions: 
<$*(r,k,,f) = ^(r,k f ,(), tf^kuqut) = <£*(r, k ; , q h t), 



^r,k ; ,ft,*), ^f(r,k,-,0 



-^7(r,k,-,f), ^(r,k,, 9i ,0 = -^(r,k i)?i ,0 and 

j£^*(r, k,, ft, i) = J^f(r, kj, ft, t) we form a coupled equa- 
tion for the quantized electric field 

J [a f (ki) + a b (ki)]^/(r, k„ t) dk, 



+ j [a°(ki)]^(r,k,-,t)dkj 



kyf— VI 



+ 2 [ J[JW*<. 9d ± ^(fc. ftW^Cr, k„ f) dkidft 
+ ff [A;. , (k,-,ft)]^°(r,k,^ / ,r)dk,d<?/ 



H.c. 



kyf— VI 



) + P(k;)]#;/(r,k;,r)dk; 



J [fi (k,)]<5?(r,kf,f)dkj 

■— ve 

JJ [fi?(k;, 9;)]^(r, ky, ft, dkydft 



and another for the quantized magnetic field 
J [a f (ki) - a\k^J\x, k„ dk,- 

+ J [a°(k i )]Jif; i (r,k i ,t)6k i 

jive 

2 [ JJ t^r (k ; , ft) T A*(ki, ft)]^ n +/ (r, k„ ft, f) dk ; dft 
JJ" [A°(k«, ft)]^, + °(r, k„ ft, f) dk,dft 



ky t - 



+ H.c. 



kyj—ve 
k r ;ive 



W(kj) - b°(kj)]Jf D Y(r, k h t) dk; 



j [b (kj)]jr p y(r,kj,t)dkj 

k x j ± ve 

£ [ jj [B f r (kj, ft) T #(k/, ft^/Cr, kj, qj , t) dkydft 

JJ [#(k;, ft)]^+°(r, ky, ft, t) dkydft 



*,,j-VC 

£ ,±ve 



Here the top (bottom) sign of the symbols ± and =p cor- 
responds to TM (TE) radiation. These two equations form 
the basis on which the coupling between the different exci- 
tations can be calculated. However, first they must be bro- 
ken down into a more convenient form. Taking Eq. ( fT5l ) and 
post-multiplying both sides by xJf T^r, k ( ', f), then integrat- 
ing over flfjf -x dfdydz (using the orthogonality relation given 
in Eq. ([8} to drop terms associated with orthogonal wavefunc- 
tions) and finally relabeling k^ — > k,, selects out the coupled 
annihilation operator equation 



ft) ± %(uj, qj)]C J r ' p (kj, qj - k,)dft, (17) 



where the coupling coefficients are given by 



Ok>ft;k;,ft) 



(18) 



A/' 



fff [%f(T,kj,qj,t) x .^(r.k.ft.f^-xdfdydz 
AT fff [E;j (r, kj, qj , t) x V7/(x, k,-, ft, f)] • x dfdydz, 



with Af = (2jT 2 hu)i)~ l . Note that we have introduced a vector 
ujj = (a>i, k y j) in Eq. ( fTTT i. This is due to the k x j component 
of the wavevector k, in the operators being transformed into 
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the frequency domain during the evaluation of the coupling 
coefficients. For SPPs in a given region, once k y is set, u> 
determines k x via the dispersion relation. For radiation the 
same is true, once k y and q are set, a> determines k x via the 
corresponding dispersion relation. Therefore k x is no longer 
a free parameter for the operators. The full analytical form 
of the couplings and further details on the domain transfer are 
provided in Appendix B. 

Next, we take Eq. ( fT5l ) and post-multiply both sides by 
Xi^j (r, k[, q[,t), integrating over JjJ -x dtdydz and finally 
relabeling — > k, and q' t — > q,. This selects out the coupled 
operator equation 

if 



[A J TM (u)i, qd + A b TM (ijJi, qi)]C'^ m (ki, qf, k,, qd 
= [b f (uj) + b\uj)]Ci i TM (k j -,k i ,q i ) 



[#(w ; , qj) ± #(w >t 9j)]C/Uk 7 -, g y ; k„ ?! )d?/19) 



Post-multiplying both sides by X J^L/(r,k[,g[,f) instead 
leads to 

[A4(M, g,) - ^(w,-, qi)]C' T ' ETE (ki, qr, k/, g,-) 
= [g/(u; J ) + S 6 (u; J )]C^(k ; ;k, >9i ) 



+ 2^ j W (Wy, ± VfrUj, qj)]C J r l m (kj, qj; k„ qddqj(20) 

Then, we take Eq. ( fT6] > and pre-multiply both sides by 
<o~f(r, k', f)x, integrate over fff xdfdydz and relabel k' — > 
ky. This selects out the coupled operator equation 

+ J] f [A{(u h qd T 9i)]C^*(ky; k„ qddqt. (21) 

Next, we take Eq. ( TTol l and pre-multiply both sides by 
S~J.(r, k'., g'., f) x , integrate over fff -x dfdydz and relabel 

*™ 7 J J 1_ JJJ 

k^ — > ky and q'j — > gy. This selects out the coupled opera- 
tor equation 

[BL(Vj, qj) ~ B b m {oJ h qj)]C J T J M * TM (kj, qf, kj, qj) 

+ J] f [#( w<1 g,) t A*( Wi> <?i)]Ok„ g/, k„ gy)dg,(22) 

Finally, pre-multiplying both sides by S~j (r, k^, q'-, t) x in- 
stead leads to 

[S4(wj, qj) + B^(u:j,qj)]C{i^(kj, qf, k h qj) 
= [af(ud-& b (ud]C J 7 i* p (k j ,qf,k i ) 

Eqs. ( fTTI i and (fT9]l-(|23T> represent six coupled equations that 
relate all the annihilation operators corresponding to the exci- 
tations involved in the scattering process. From these equa- 
tions we are able to form a transfer matrix, as depicted in 



Fig. |2](a). This matrix will then provide a complete descrip- 
tion of the quantum dynamics for the scattering process that 
occurs at the interface. Note that on the right hand side of 
each of the six coupled equations we have kept coupling co- 
efficients which, by definition of Eqs. © and ( fT8l ), should 
strictly be equal to unity. We keep these coefficients in order 
to allow for the inclusion of evanescent excitations in our the- 
ory which have a normalization equivalent to Eq. ([8]l up to an 
overall factor ±i 13511 . which we will specify when required. 



V. TRANSFER MATRIX 

In order to form the transfer matrix we take Eqs. dTTb and 
([T9l>-(l2"3ll and use truncated summations to approximate the q t 
and qj integrals for the radiation excitations H 1 3fl - For accu- 
racy and convergence Gaussian quadrature summation [36] is 
employed, where an arbitrary integral can be written as 



b b 

f(q) dq = — 



f 1 b - a N+l 

■ J f(u) du = — — V w m f(u m ) 



(24) 

with u = 2q/(b — a) + (a + b)/(a — b), the weights w m - 
2/[(l — u^j){P' N +l (u m )) 2 ] and the u m abscissa chosen to be 
the roots (zeros) of the P^+\{u) the Legendre polynomial. In 
our case we have the integral limits a — and b = q cut . 
Here, care must be taken with the weights. We must en- 
sure that the sum of the modulus squared of a given set of 
amplitudes, for example a£(u>;, u m ) (= a£ (a?,-, qf)), calcu- 
lated from the coupled equations, for the discretized radia- 
tion operators A*(a>y, u m ) (= Aj! m (w;, g,)) leads approximately 
to the total back-scattered energy, i.e. J \a h r (a;,-, g/)| 2 dg,- ~ 
2m=i l«r,m( w i.9i)| 2 - But > following Eq. ([24), we have the re- 
lation J \c4(Wi,qd\ 2 dqi = Y!Lt\ w ' m l^mC^i.^i)! 2 , with w' m = 
Wm gcut/2. Thus we rescale all operators involved in integra- 
tions, e.g. A b r (u„ qi) — > A b rm (u>i, qrj *jw^ l to compensate and 
write the six coupled operator equations in discretized form 
as 

[af{Ui) + a V0]C/k<; ki) = [b f {ujj) + b b (lj ^{kf, k,X25) 



Ar+i 



+ J] iHjUj, qj) ± g/)K',';,(k, u m ; k,) J<, 

r,m= 1 

[A f TK „(u>i, qi ) + A b M n (u)i, qi)]C'^ m (ki, u n ; k,-, u n ) = (26) 
[BfQaj) + S V;)]C ; f™(k/; k ; , u n ) V< 

w «' ^i' ^n) 

r 

N+l 

+ $ f r,m(Uj r qj) ± B b } . m (uJ h q j)]C J r l m (kj , U m \ k/, U n ) ^W f m W f n , 
r,m=\ 

[A{ Kn (Ui, qd - A b Tm (u)i, qi)]C^ m (ki, u„; k ; , u n ) = (27) 

[Bf(wj) + &(.uj)}P*J)Lfi k h u n ) V< 
+ V [B{„(LJj, qj ) ± B b n (ujj, qjjlC^ikj, u n ; k h u n ) 
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+ Yj ${m(Uj,qj) ± B b rm (u)j,qj)]C J r ' TE (kj, u m ; k„ u„) y/w' m w' n , 



r,m= 1 
Iff. . , tb 



W{uy)-b°{uyj]C pp (k/,k 7 ) = 



(28) 



i 

+ [MA^i^i)TA h rm {^i,q i )]Cj,' r \kj\k i ,uJ^J^, 



N+l 



r,m= 1 
if 



[Bmjuj, qj) - B b TKn {u> h qj)}Ci ] M * TM (kj, u n ; k,, u n ) = (29) 

[^(w,) - aPiudV&QLj, u n ; k,) V< 
+ ^}Kn(Ui, qd T ^(u;,-, ^)]C^*(k y , m„; k„ «„) 




b) 









b 6 



AT+] 



+ 2 t^{m( a; ''' id T A h rm (u)i, qi)]C J T l*.(kj, u„; k„ « m ) V<X> 

r,m— 1 

[#i,n(^./, ?;) + B b TE „{u: h qj)]C J T ] E * E (kj, u n ; k jt u„) = (30) 

-ta^u;,) - a*(w,)]C^;(k Jt «„; k,) V< 
+ ^J.M,n(Ui, qd T A b r „(u) b qi)]C J TE *(kj, u n ; k h u„) 

r 

N+l 

+ Yj &JM> 9d T ArJ-Vt' qi)]C ] r E *(k h u n ; k„ u m ) y/w' m w' n , in Fig.|2](a), for the operators at the interface as 



FIG. 2: (Color online) Transfer matrix and SPP beamsplitter angles, 
(a): Input-output transfer matrix T linking all the excitations at the 
interface, (b): Angles involved in the scattering process: the incident 
angle 9^ (#,■ ) of the SPP excitation a 1 (b b ) in region i (J), the trans- 
mitted angle (?,-, (6j t ) of the SPP excitation V (a ) in region j (;). The 
relations (?,,. = 0j, and 6j r = must hold for the excited outputs to be 
indistinguishable. This can always be achieved, see text for details. 



r,m= 1 

where we have chosen u m (u„) to represent qj (qf) on the first 
three equations and vice-versa for the second three equations 
in the coupling coefficients. Doing this allows us to write 
the coupled equations in a more compact form as four matrix 
equations 



( A b \ 

TM 


1 


B f 




A b 




•. B TE , 


\ 



T31 



T22 
T32 
T42 



T13 

T23 
T33 

T43 



T 14 ^ 
T 2 4 
T34 
T44 



B' 



(36) 



where A f r ' b := (A™, 



if* 



Here 



N+l 



& 4- A b 



r m=0 
N+l 



ftf ftb 

D m,n D m,n 



&f — A b 

TEJl TE,n 



r m=0 
N+l 



J]Y(B{ jn ±B b .J(D rTE )„ 

r m=0 
N+l 

Bte,h + B TEn — *y ^ *y ^ (A rm A rm )(F TEr ) m 



(31) 
(32) 
(33) 
(34) 



r m=0 



Here, Eqs. (|25) and (|26) lead to Eq. CD}, Eqs. ([28} and (g9) 
lead to Eq. ( T32l , Eq. ( l271 i gives Eq. ( l33l l and Eq. ( f30T > gives 
Eq. (O, where 



(35) 



i^„:={a / -Vol« = o}, 

A^:={0|« = 0}, 

:= {^(w,, ft) = 4'*(w ;> «„)| n > 0}, 

^:={0|« = 0}, 

^ := {B{„ & (^, ?y) = u n )\ n > 0}, 

with the D and F matrices given in Appendix C. From 
Eqs. (l3Tb-(f34T> we can finally write the transfer matrix, shown 



) and similarly for B r ' . 
pj ,,, At n , „ and B b „ represent the forward and back- 

™,0' ™,0' TM,0 TM,0 r 

ward propagating SPP excitation operators, as specified in 
Eq. d3~5l l. To obtain the full transfer matrix T, one must calcu- 
late the elements of the D and F matrices, which consist of the 
various coupling coefficients C^ v (kj, qf k,-, qi) (see Appendix 
B). The matrices D and F then enter into the elements of T as 
outlined in Appendix D. 

One of the fundamental differences between the scattering 
of SPPs described here and the scattering of photons at a stan- 
dard optical interface is that the spatial profile of an incoming 
SPP's wavefunction is modified due to the variation in the per- 
mitivitty of the materials on either side of the interface (metal 
and dielectric). As a result, the reflection and refraction of 
SPPs do not follow the standard laws of optics for determin- 
ing the various scattering angles involved, such as Snell's law 
and the usual Fresnel equations 1IT4I fl9[l . Fortunately, all the 
information required to obtain the scattering angles for the 
SPPs and radiation excitations is contained within the cou- 
pling elements. In particular, one can immediately and in a 
straightforward manner calculate the angles for all the SPPs 
involved in the scattering process. 

As can be seen from Eq. ( IB- II) for the coupling coefficients 



Cjl v (kj, qy, k,-, q/), the magnitude of the wavevector in the y di- 
rection, k y , and frequency of the excitation, a>, must be the 
same on either side of the interface. This allows one to im- 
mediately specify reflected and transmitted angles for all SPP 
excitations involved. This can be achieved because for SPPs, 
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once a> is set, the DR given by k = (a)/c)(^e m /(Q + e m )) ' au- 
tomatically sets the value of k, and with k y also set one directly 
obtains the value of k x . For instance, consider a single-SPP 
excitation in region ; incoming at specific angle 8/ i with a set 
frequency <j;. Recently it has been shown that such a quantum 
excitation can be efficiently generated on the metal-dielectric 
interface using an attenuated-reflection geometry [7, 8]. With 
all the material properties of the interface configuration set, 
using the DR for SPPs in region i, the frequency and angle val- 
ues result in a corresponding wavenumber k t and wavevector 
k, = k t i x + ky$. From this, we can calculate the wavevector 
for the transmitted SPP in region j, as the condition loj = <w, 
and k y j = k y j must be met across the interface. From the 
wavevector, the transmitted angle 0,, can be found. Note that 
due to the monoticity of the DR, there is only one SPP in re- 
gion j that can be excited by a single SPP incident in region i. 
Furthermore, based on the above arguments, it is straightfor- 
ward to check that the incident angle is equal to the reflected 
angle, i.e. Q k = Q lr . For the radiation excitations such a direct 
specification of the angles, as in the SPP case, is not possible. 
This is because k (and therefore k x ) depends also on the value 
of the parameter q. Thus, for a set a> and k y there are a con- 
tinuum of radiation excitations which couple to an incoming 
SPP, each one exiting the scattering region at a particular an- 
gle with a given probability amplitude defined by the transfer 
matrix T. 



VI. BEAMSPLITTER CHARACTERIZATION 

A beamsplitter is one of the most important devices in linear 
optics. It constitutes an essential building block in any classi- 
cal optical setup [37] and in the context of QIP it plays a cru- 
cial role in optical implementations of quantum gates ll38i[39ll . 
In the simplest case of a lossless optical device, the action of 
a photonic beamsplitter can be described mathematically by a 
2x2 matrix that relates the two output ports to the two input 
ports. It is possible to derive further constraints and relations 
for the elements of this matrix, the reflection and transmission 
coefficients, based on simple energy conservation arguments 
or on the fact that the device is passive, reciprocal, and lossless 
MM- 

The situation is somewhat more involved for SPPs because 
the incoming SPP excitations are not only coupled to the out- 
going SPP ones, but also to a continuum of radiation excita- 
tions. In this sense, the problem is similar to the case of a 
lossy beamsplitter and, with some care taken - using the the- 
ory derived in the previous section - it can be treated as such. 
Furthermore, we note that the geometry considered here con- 
sists of two different media on either side of the interface, 
thus unlike the standard photonic beamsplitter scenario, one 
must also take into account the different group velocities as- 
sociated with the propagating plasmonic excitations. For the 
present discussion we will be interested in the main beam- 
splitting properties in the immediate vicinity of the interface. 
Propagation issues will be discussed in the next section. 

It is also important to note that while a beamsplitter is con- 
ventionally defined as a device used to divide a single incident 



beam (on one side) into two beams (one on either side) with 
a ratio of intensities, in this work we focus on the description 
of a beamsplitting device where SPPs can be incident on the 
interface from one side or the other, and even from both sides 
at the same time, but with the condition that the output modes 
for the excitations are the same, regardless of which side the 
SPPs are incident on. This scenario is the usual symmetric 
case considered for a beamsplitter operating at the quantum 
level which mixes two input fields and produces two output 
fields lf20L l40h - In order to highlight that the device must be- 
have in this symmetric way we focus our discussion on the 
latter case - when a single-SPP is incident on both sides. One 
can easily modify this for other cases, such as when a single- 
SPP is incident from only one side or the other by placing one 
of the SPP excitation field modes in the vacuum state [20]. 

First, consider that an SPP is incoming from region j at 
an angle 0,- and another is incoming from region i at an an- 
gle 6i r In the present context, for a quantum beamsplit- 
ter operating faithfully at the single-SPP level it is essential 
that the input SPP excitations produce outputs which are in- 
distinguishable in both their spatial and temporal degrees of 
freedom, regardless of the transmission and reflection coef- 
ficients. One major requirement for spatial indistinguisha- 
bility is that the angle 6 Jt , corresponding to the transmitted 
SPP into region i, must be equal to the angle 6 h of the re- 
flected SPP excitation incoming from region i. In addition, 
we must ensure that 0j r = These conditions are shown 
more clearly in Fig. [5] (b). Surprisingly, such a matching 
set of splitting angles can be achieved. By setting the inci- 
dence angle for the incoming SPP in region j to 6j t = Q h - 
cos" 1 [(1 - e d je,„j(e dJ + e mJ ) sin 2 U I 'e^fe,; + e m>i )e m j) 1/2 ], we 
find that the transmitted angle 6j, = Q h regardless of the fre- 
quency at of the excitation. Thus one is always able to achieve 
the required matching for the angles in the scattering process, 
as long as 6^ results in g [0,7r/2], otherwise plasmonic to- 
tal internal reflection (TIR) will occur. This symmetric split- 
ting regime has not been investigated before, even in the clas- 
sical domain. In what follows in our beamsplitter characteri- 
zation it should always be assumed that the matching angle set 
has been chosen as specified above, regardless of the param- 
eter regime being studied. All angles will be specified where 
relevant. 

For temporal indistinguishability, it is important that the in- 
coming SPPs have exactly the same arrival time at the inter- 
face, as well as having the same spectral profile (bandwidth). 
For the former, we assume this can always be achieved ex- 
perimentally, constituting a practical issue rather than any- 
thing fundamental in the theory. For the latter, we will begin 
our discussion with a simple single-mode picture and discuss 
wavepacket excitations later. 

Continuing with our aim to make the scattering interface 
work as an SPP beamsplitter, we must also quantify how well 
SPP excitations couple to each other and not to the radiation 
excitations. From the transfer matrix given in Eq. fl36l l one 
can link the output to the input operators for the SPPs using 
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the relations 



N+1 



(T 11 ) 00 a / (w ( ) + YjFidQeALjPhqd (37) 
e=i 

N+l 

+(r 12 ) 00 &(uj) + Yj(T lz ) oe B b me (u>j, qj) 

e=\ 

N+l N+l 

e=i (=i 

N+l 

(r 21 )„ () + 2( r ai (38) 

e=i 

N+l 

ot^TMii^p qj) 

1=1 

N+l N+l 

+ 2(r 23 )oM4 / (o; / , q,) + YJJtaKx&te^p 9j)- 
t=i t=i 

Similar relations can be written for the radiation excitations. 

In order to achieve perfect matching of the quantized elec- 
tromagnetic field at the interface, we have included evanescent 
radiation excitations in our theory and therefore increased q cut 
from tOi(edj) 1 ^ 2 /c to a>i(edj — e m j)^ 2 /c. These excitations rep- 
resent a contribution to the total field which decays exponen- 
tially from the interface in the x-direction (in addition to the 
usual decay in the z-direction parameterized by the variable 
q, as described in section II A). Although these excitations 
do not propagate they are important in the study of near-field 
phenomena. For instance, such near-field radiation could be 
exploited to achieve coupling with a resonant auxiliary sys- 
tem placed in close proximity to the interface. However, by 
extending the quantization procedure followed for the propa- 
gating radiation (see Section II) to decaying evanescent radia- 
tion, we have that the photon number operator of the evanes- 
cent excitations is not Hermitian. This implies that such an 
operator is not a physical observable. Nonetheless, when we 
look at the scattering problem within the far-field approxima- 
tion we have that the transfer matrix T is indeed unitary and 
provides amplitudes for the transfer of the physically well- 
defined propagating excitations. A more rigorous treatment of 
the evanescent excitations along the lines of Ref. [411 is pos- 
sible, but such a treatment is beyond of the scope of this work. 
Moreover, the appearance of evanescent excitations and their 
non-reciprocal behavior on either side of the interface leads to 
a general lack of reciprocity in the splitting device, the effects 
of which can be measured, as described next. 

From Eqs. ( f37T > and fl38T > we define the coefficients for the 
SPP beamsplitter of an incoming SPP in region i as follows, 
transmission r = |(7"2i)oo| 2 , reflection p = \(Th)qq\ 2 , and to- 
tal loss (into propagating radiation) <x = |(rii) ro | 2 + 

„mmaxj+l i,_ . i 2 ^mmaxi+1 \/r r \ \2 , ^mmaxi+1 \ frr \ \2 

Li l= \ |U2i)ffl| +Le=i |U3i)fo| +L( =J ' \{i4i)m\ , 
where mmaxj {mmaxj) stands for the number of non- 
evanescent (propagating) radiation excitations in region i (j) 
upon discretization. Similarly we define the coefficients of an 
incoming SPP in region j as t' = |(7"2i )oo | 2 for transmission, 
pf = |(7ii)oo| 2 for reflection, and a' = ^'"T^ \( T n)eo\ 2 + 




0.6 0.9 1.2 

6 H (rad) 



FIG. 3: (Color online) Scattering interface as an SPP beamsplitter 
for /lo = 790nm shown in panels (a) & (b) and for 1500nm in pan- 
els (c) & (d). (a) & (c): Forward SPP transmission r (solid lines), 
reflection p (dashed lines) and scattering cr (dotted lines) coefficients 
versus the incidence angle In all plots ed,i increases in steps of 
0.25 from 1 to 3 going from right to left (bottom to top for the scat- 
tering coefficient), (b) & (d): Coefficients |t — t'|, \p— p'\, \cr — cr'\. 
In all cases, e dJ = 1 and a pM = cof rer = 1.402 x 10 l6 rad/s are 
chosen. Note that in the plots of panels (a) & (c) the coefficients r, 
p and cr are limited at large angles of d\ i . This is due to a reduction 
in numerical precision close to plasmon TIR for which we do not 
include the corresponding values. 



x-\tnmaXi+\ \/r r s i 2 ^mmaxj+l \ /rr * n ^mmaxj+i \,r r 
hf=\ \K 1 2Vm\ + 2jf = i |U3l)fl)| + Zif=l K-Mjflll 

for total loss, which are calculated for the reverse configura- 
tion (from j to i). Thus, using the above relations we can 
rearrange the transfer matrix Eq. d36*l > as 



' a b (uii) ) 


i 






■. '■ 


V 



a f (u>i) 



, (39) 



where the phases of the coefficients stem automatically from 
the matching of the transverse components of the wavefunc- 
tions in Eqs. $13[ and (TBI) . We find that when e^j > e^j me 
phase <p — and when ejj < €dj the phase <p — n. In order to 
check numerically how close the interface devised here is to a 
reciprocal one, we require the conditions t ~ t', p ~ p' and 
cr ~ cr' to be met H20il . 

In Fig. [3] we show how the SPP transmission, reflection and 
scattering coefficients vary as a function of the incidence angle 
Q u and the value of e^j for free-space wavelengths Aq = 790nm 
(a) & (b) and 1500nm (c) & (d) (w, = 2nc/A ). These wave- 
lengths have been chosen as examples to illustrate our results 
as they correspond to those used regularly in quantum op- 
tics experiments 13711 . The range of €dj : 1 — > 3 (in steps 
of 0.25 going from right to left) chosen corresponds to that 
which can be obtained with the use of basic optical mate- 
rials ll42ll . Here, as an example we set e^j = F with the 
metal in both regions modeled as silver. We also choose 
the number of radiation excitations to be N + 1 = 200 and 
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q cat = 10max[(Q,^ - kj/^Ae^kj - ^'jHn order to 
satisfy convergent behavior of the coefficients 01311 . One can 
see from panels (a) and (c) that, in the case of normal inci- 
dence 6j j = 0, a large value of e^j — 3 is required in order to 
observe a significant departure of the transmission coefficient 
r from unity. However, this is accompanied by a significant 
increase in the value of cr corresponding to scattering into ra- 
diation excitations. To avoid such a scenario one could use 
lower ed/s and vary the angle 6*, ; instead. However, one must 
be careful in doing so for two reasons. First, there exists a 
critical value of the incidence angle for which plasmonic TIR 
takes place. In order to ensure that we do not reach the plas- 
mon critical angle our computation stops at those values of 
Q Vi near which TIR for the radiation excitations occurs. This 
can be seen more clearly in Fig. [4] where we show the depen- 
dence of the transmission angle (as it approaches n/2 rad) 
as a function of the incident angle 6 ti and e^,-. However, it is 
the stable region around this critical angle that there is a wide 
range of opportunity in setting the transmission and reflection 
coefficients. The second reason one needs to be careful is di- 
rectly related to the first, and it is that there is a rapid change 
in the transmission and reflection coefficients in the region of 
opportunity, as shown in Fig. [3] (a) and (c). Thus, the coeffi- 
cients are quite sensitive to the incidence angle #, f . In addition, 
Fig. [3](b) and (d) one can see that the difference between the 
coefficients, i.e. |r — r'|, \p — p'\, \cr — cr'|, is relatively small 
with varying 6j j and edj- 

In our calculations we have set e c /j = 1 for simplicity. In 
general, one needs to be careful with the opposite configura- 
tion, where the incoming excitation goes from a higher per- 
mittivity region to a lower one. This is because the coupling 
to evanescent radiation appears to be more substantial. We 
also note that we have verified numerically that when metals 
with highly negative permittivity are used (challenging from a 
fabrication viewpoint j43ll44ll ). compared to the values of the 
dielectrics, it is possible to favor the stability of the SPP ex- 
citations over the radiation, suppressing losses into radiation 
even more. 

While we have considered the dielectric media and metal 
in both regions to be passive (fixed) in this work, our theory 
applies equally well to the case of active materials, such as 
electro-optical polymers or voltage programmable liquid ll45ll 
for the dielectric media, and metal-semiconductor quantum 



(a) (b) 




B H (rad) 9 n (rad) 



FIG. 4: (Color online) Transmission angle 8j, as a function of the 
incidence angle 9t t . (a): Ao = 790nm. (b): Ao = 1500nm. In both 
plots £d,i increases in steps of 0.25 from 1 to 3 going from right to 
left. 



well structures [46] in place of the metal. The use of such 
active material would provide greater flexibility and external 
control of the parameter regime of the constructed beamsplit- 
ter, and thus the range of transmission and reflection coeffi- 
cients, without the need for re-fabrication of the device. More- 
over, such 'on-chip' beamsplitters may be created and de- 
stroyed at desired locations using the correct associated elec- 
trical circuitry lfl4ll . 

In summary, we have found that operating close to the crit- 
ical angle for 6^ provides great flexibility in obtaining a range 
of values of the transmission and reflection coefficients, how- 
ever it comes at the expense of the reciprocity of the beam- 
splitter, loss into radiation, as well as the necessity to accu- 
rately set the angle #, ; due to the rapidly changing values of 
the coefficients. As is evident from the above initial investiga- 
tion, the physical properties of the materials involved together 
with the incidence angle determine how well one can optimize 
an SPP beamsplitter to provide a desired splitting ratio. 

While we have found that a range of splitting ratios can 
be reached, as an example we now fix the goal of achieving a 
50:50 plasmonic beamsplitter, a versatile component in classi- 
cal and quantum optics 13711 . By allowing the different param- 
eters of the interface to vary, such as the dielectric material, 
metal, frequency and angle of incidence, we seek to optimize 
the beamsplitter's performance. 

In Fig. [5] we show the numerical results for an optimization 
procedure we have followed in order to bring the scattering 
configuration as close as possible to a 50:50 beamsplitter for 

= 790nm (a-f) and 1500nm (g-1). The intent of these plots 
is to show that as the angle 6 it is varied, the dielectric con- 
stants and plasma frequencies of the metals on either side can 
be changed in order to bring the beamsplitter as close as pos- 
sible to a 50:50 splitting ratio. Here we have imposed several 
constraints on the acceptable configuration. First, the maxi- 
mum amount of scattering into radiation cr should be 5% or 
less and second, in order to force the interface to work as a re- 
ciprocal device, we constrain the differences |r — t'|, \p — p'\ 
and \cr — a J \ to be lower than 2.5%. One can clearly see in 
Fig. 0(a) ((g)) that a 50:50 beamsplitter can be reached for 
Aq = 790nm (Ao = 1500nm) if we allow the metals to vary 
across the interface as shown in panel (b) ((h)), assuming 
u>p lver < co p j/j < 2oj p dver , for a range of incidence angles 6 ij . 
Such a range in the value of the plasma frequency for the met- 
als should be possible, for instance by embedding highly con- 
ducting surfaces perforated by holes l43tl or depositing pat- 
terned metallic films on semiconductor structures 14411 . On 
the other hand, if we set the metal as silver on both sides of 
the interface as a less demanding scenario as shown in panel 
(e) and (k) for Ao = 790nm and Ao = 1500nm respectively, 
we can reach a 50:50 beamsplitter, although only at the higher 
wavelength of Ao = 1500nm, as highlighted by panels (d) and 
CD- 

Throughout the above analysis we have assumed that the 
scattering of SPP excitations into radiation excitations can be 
treated as a loss mechanism. In order to ensure that this is in- 
deed the case and that the radiation excitations are not able to 
interfere with any SPP detection process or subsequent SPP 
beamsplitter operations one must check that the radiation ex- 
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FIG. 5: (Color online) Results of the optimization procedure for achieving a 50:50 SPP beamsplitter with the scattering interface over a range of 
incidence angles ft for Ao = 790nm (1500nm). (a) ((g)): SPP transmission t (solid), reflection p (dash), and scattering a (dotted) coefficients, 
(b) ((h)): Values of e d ,i (solid), ft, (dash), (j Pii (dotted), w pJ (dash-dot), (c) ((i)): Coefficients |t — t'| (solid), |p -p'| (dash), \cr - tr'| (dotted), 
(d-f) ((j-l)): Same as (a-c) ((g-i)) but with the restriction a> p j = co p j = 1. In all cases, ejj = 1 is chosen and values of a> p j/j are measured in 
units of ai'p lve ' = 1.402 x 10 16 rad/s. The vertical dotted lines designate optimized angles and the corresponding material parameters for the 
50:50 beamsplitter. 



citations scatter in a direction such that their quantum degrees 
of freedom can effectively be discarded, or 'traced out' 1381 . 
from the SPP dynamics. In Fig. [6] we show the scattering 
angles of the resulting TM radiation excitations from an in- 
coming SPP excitation in region i with A$ = 790nm ((a)) and 
1500nm ((b) and (c)). The parameters chosen for the mate- 
rials and the incidence angles correspond to those defined in 
the caption by the respective dotted lines of Fig.|5]and given in 
the table of Fig. [6] For these three scenarios, the scattered TE 
radiation is an order of magnitude less than the TM radiation. 
As the total radiated power has a maximum of 5%, we focus 
on the TM excitations as the main sources of power loss. The 
plots of Fig. [6] essentially show the direction of the 2(N + 1) 
scattered radiation excitations - defined by the wavevectors 
(k X) i, kyj, qf) and similarly for region j - multiplied by their 
contribution to the total scattered power (normalized by the 
largest contributor). The start (end) points of these normal- 
ized 'power- vectors' lie at the origin (on the curves). For both 
= 790nm and 1500nm one can clearly see that the result- 
ing radiation excited in both forward and backward directions 
exits the scattering region at reasonably large enough angles 
from the x-y plane such that it is possible to trace them out 
from the system dynamics and thus treat the radiation as a 
truly lossy mechanism. Similar plots can be made for the SPP 
incoming from region j with the same parameters and result- 
ing in the same conclusions. Note the fact that we can trace- 
out the radiation excitations is not a general rule even though 
it applies to the parameter range optimized and investigated 



here. One should analyse the TM and TE radiation scattering 
on a case-by-case basis. 

VII. QUANTUM INTERFERENCE 

In our derivation of the properties of the scattering inter- 
face, we have worked explicitly with the quantized form of 
SPP and photon radiation fields. We have done this because 
at the quantum level there are interactions between SPP states 
which cannot be described within a classical framework. In 
this context a particularly enlightening scenario is the inter- 
ference of two single-SPPs. We now apply our theory to this 
case as an explicit example of the necessity of a quantum the- 
ory for plasmonic scattering. 

Consider an input state with two SPPs, one on each side of 
the interface. For simplicity, we assume that both SPPs are 
monochromatic, with the same frequency, and arrive in co- 
incidence. Thus temporal indistinguishability is immediately 
satisfied. Although, more realistic Gaussian wavepackets will 
be considered in a subsequent study with more complex scat- 
tering geometries 04711 . in the present work we have verified 
numerically that when a lOnm- width wavepacket centered at 
the frequencies corresponding to Aq = 790nm and 1500nm 
are considered, the values of the coefficients of the scatter- 
ing matrix do not change appreciably. Thus a single-mode 
picture can be adopted as a close approximation for this ini- 
tial investigation. We also note that the spectral shape of a 
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wavepacket does not change on passing from one region to 
another. The situation changes during propagation, where af- 
ter a given distance, the group velocity related to two different 
media causes a deformation (broadening) of the wavepacket 
at different rates [37]. This can be compensated with the in- 
corporation of appropriate dielectic media on the metal sur- 
face ll47ll or by considering only short propagation distances. 

We write |1,0) 1 1, 0)^, . r . in order to represent the input 
state in the Schrodinger picture, which consists of a single 
SPP on each side of the interface plus vacuum for the input 
radiation. To calculate the output after scattering occurs at the 
interface we apply the transfer matrix in Eq. (1391 to the input 
state, written in the Heisenberg picture as 

1 1, o) Piri \i,o) Pjr . = aVfaWfrj) |o, o) Piri |o, o) m . (40) 

Due to the non-zero value of the couplings between SPP and 
photon radiation, the output state will contain terms in which 
the radiation excitations are populated with a given probability 
amplitude. More precisely, for the input defined in Eq. d40b we 
have the output 




FIG. 6: (Color online) Scattered TM radiation plots resulting from 
the scattering of an incoming SPP from region ('. (a): 50:50 beam- 
splitter for Ao = 790nm. Here, the parameters chosen correspond to 
those defined by the dotted line of Fig.[5](a), where the metal dielec- 
tric is allowed to vary, (b): 50:50 beamsplitter for Ao = 1500nm. 
Here, the parameters chosen correspond to those defined by the dot- 
ted line of Fig. |5](g). (c): 50:50 beamsplitter for Ao = 1500nm. Here, 
the parameters chosen correspond to those defined by the dotted line 
of Fig.|5](j), where the metal is the same on both sides of the inter- 
face. As one cannot reach 50:50 for Ao = 790nm in this regime, such 
a plot has not been included. 



( V,0 - «T* y/p&'iuj) + ...) |0, 0) m |0, 0) p . r . 
= ^r{f \2,0) m |0,0) ;w - e-» \Q,0) m \2,0) p J 

+ (T-p)|l,0) wl |l,0> w + ..., (41) 

where the dots in the last line account only for terms in the su- 
perposition with strictly zero SPP excitations. Assuming that 
the state at the output SPP ports after the scattering process 
can be measured ideally i48tl . then one can calculate the prob- 
ability to measure a single SPP at each output port in 'coinci- 
dence', i.e. a joint detection event. This probability is found 
to be P coincidence = (j — p) 2 JMl- In particular, for a 50:50 
beamsplitter, a zero probability is obtained. This is a result of 
the well-known Hong-Ou-Mandel effect [49] and it is useful 
in characterizing how well quantum interference can occur for 
a particular scattering device. This effect is purely quantum- 
mechanical and cannot be described in terms of a classical 
treatment of the SPP and radiation fields H49fl . Indeed, us- 
ing a classical description of the fields leads to a minimum 
probability of 0.5 ll50ll . An experimental test of this scattering 
scenario would yield strong evidence supporting the bosonic 
nature of SPPs in the frequency regime considered. We note 
that previous experiments J3, [j, [Toil have tested some of the 
basic quantum properties of SPPs indirectly. The approach 
presented here would be more direct and conclusive. 

In Fig. [7] we show the values of the coincidence probability 
obtained after the optimization procedure carried out accord- 
ing to Fig. [5] One can clearly see how the coincidence drops 
to zero for values of the incidence angle approaching a 50:50 
beamsplitter. 



VIII. SUMMARY 

In this work we introduced the quantum mechanical formal- 
ism for treating SPP interactions via scattering at an interface. 
To do this we extended a recent classical study to the quan- 
tum domain and at the same time generalized it to allow for 
the consideration of SPPs scattering at an angle. This work 
provides a first step in the direction of showing SPP quan- 
tum interference in addition to understanding SPP circuitry 
at the quantum level. Here, we used our theory to investi- 
gate the possibility of achieving plasmonic beamsplitters that 
operate faithfully at the single SPP-level. We found that a 
wide-range of splitting ratios can be reached depending on 
the frequency and angle of the input SPP excitations and that 
power loss due to unavoidable scattering into photon radiation 
can be suppressed to 5% or even less in some cases. As the 
beamsplitter geometries we have investigated are readily ac- 
cessible to experiments due to the basic properties of the ma- 
terials involved, the theoretical findings reported here could 
be tested within the near-future. As an application of our the- 
ory, and an example of one such possible test, we investigated 
theoretically first-order quantum interference effects of SPPs 
using a 50:50 plasmonic beamsplitter optimized and incorpo- 
rated into a Hong-Ou-Mandel type setup. Our study is the first 
to show that surface plasmon beamsplitters can reliably oper- 
ate at the quantum level. Furthermore, our developed theory 
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FIG. 7: (Color online) First-order quantum interference of SPPs where the parameters chosen correspond to the optimization procedure used 
for the 50:50 SPP beamsplitter described in Fig.[5]for Aq = 790nm (1500nm). (a) ((c)): Two-fold coincidence probabilities, (b) ((d)): Two- 
fold coincidence probabilities with the restriction u) PI = oj Pj = aif ver . In all cases, ej ] = 1 is chosen, (e): Optical beamsplitter analogy for 
the plasmonic beamsplitter with two detectors measuring the presence of an SPP at both ports at the same time (C represents a coincidence 
detection circuit). 



could be applied equally well to other waveguide geometries 
such as long-range, channel or dielectrically loaded B16l4l8ll . 
Thus we expect our work to help open up new directions of 
research into the design of efficient and practical components 
for on-chip plasmonic-based QIP at the nanoscale. 
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APPENDIX A 

The positive and negative frequency parts of the electric and 
magnetic quantized radiation fields are given by 



E+ (k h q u r, t) ■■ 
E-(ki,qi,r,t) ■■ 
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lP r (k i ,q i ,z)A r (k i ,q i )e- i ^ t -'' > 
^:(k i , qi ,z)Al(k i , qi )e i ^ r ''> 



(cos qiz + n. sin #;z)#(z) 









T 7 ) 



and 

ip TE (k h q u z) = my? 



(iki + ^-z)e v ' z §(-z) + 

Vi 



[/k,(cos qi z — — sin qi z) + — z(cos q t z + — sin g,z)]#(z) 
Vi Vi qi 

These can be found from the relation ip r (kj, qi, z)e~ ix ^ ,,g,,T ' f) = 
V x <p r (ki, qi, z ) e *cr-to,t and^kj, q u r, t) = cojt - k, r n/2. 
For the SPP field we have 
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which can be found from the relation i/> p (k;, z)e '* (k » r >'' = V x 
4> p (ki, z ) e ik r'-^> and ^(k,-, r, t) = w,f - k, r - n/2. 

The fields for an excitation defined by k, and q t are 
given by E^k,-, q u r, t) = E" (k ; , q u r, f) + E+.(k ; , q h r, t) 

and %i(ki,qi,r,t) = ft-(k,-,g,-,r,0 + ft^(k,-, q u r, t). 
The total positive/negative fields are given by 
(2n)~ 1 jg* dkid qi E±(ki, q h r, t) and 

= (2n)- 1 ffi™ dk^H^Ck,-, gi , r, f). The 
total fields are E^,(r, t) = E+(r, t) + E~.(r, t) and 
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APPENDIX B 

For convenience we provide a compact analytical form for 
the coupling of Eq. ( fl~8T > given by 
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Here we have r mi (k h qd = (iedjVj - e^qd/fieajVi + e mii qd, 
as well as r TEi (kj,qd = (zV; — qd/(ivi + qd, with similar 
definitions for side j together with the function 6+ (k) = 

I- h{z)e' kz dz = \5(k) — ^yy, where /z(z) represents the 
Heaviside function, V stands for the Cauchy principal value 
of the integral and 6(k) = ^ J ^ e' kz dz. 

Note that each of the above couplings will be integrated 
over k x , k y and q for the relevant region due to the integrals 
appearing in Eqs. (fl~5T > and ( TToT l. This gives physical mean- 
ing to the delta functions which select out the operators spec- 
ified in the coupled operator equations, Eqs. (fTTT i and (fT~9T >- 
(l23l l. Here, the integration over k x is changed to integration 
over a) for convenience. For example, for SPPs we have 
dkxj — > (ki/k X jVc,ddu>i with the operators transforming ac- 
cordingly, aHkd — > (k x jVcj/kd l ^ 2 af(u>j,kyj). Thus, all oper- 
ators are strictly transformed from wavenumber to frequency 
domain for the k x component of the wavevector, hence the no- 
tation af{ujd f° r the operators, where u;; = (tOi,k y j). 



APPENDIX C 

The D and F matrices appearing in the set of coupled ma- 
trix equations (f3TT)-(f34b are defined in terms of the elements 
of C matrices, which are themselves defined in terms of the 
couplings Cu V (]£j, qf, k,, qd (specified by Eq. ( TT~8T > and given in 
Appendix B) as follows: 



(C m 7 M )oo — Cpp(kj> kj), 
(C T MTM)on — Cp m (kj', kj, u„) yjw' n , 

(C I M I 'M )m() = Cmp(kj, u m ; kj) ^Jw^, 

(C TMTE )on = C J p ' TE (kf k it u„) V<, 

{C r r)nn ~ Cl r (kj, U n , k/, Ufd, 

{C rr )mn = C£(k j, U m \ k/, U n ) yjw' m w' n , m + U, 



(C-l) 
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with zero elements 

(C TETM )mn = Vm,n, 

(CVe7"e)()0 — (C^^O/i — (C TETE ) m Q — 0. 

From the above matrices we then have the D matrices 
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Here mmaxt (mmaxj) represents the number of non- 
evanescent (propagating) radiation excitations in region i (j) 
upon discretization. The imaginary factors account for the dif- 
ferent normalization of non-evanescent and evanescent exci- 
tations, as described in the main text. Note that the subscripts 
TM and TE denote a different D or F matrix, whereas the in- 
dices m and n denote the elements of these matrices. In the 
limit of normal incidence, i.e. 6, j = 0, coupling to TE ex- 
citations does not occur and the D and F matrices lead to a 
transfer matrix as given previously in Ref . lfl3ll . 

APPENDIX D 

The T/j entries of the transfer matrix are given by 
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and the F matrices 
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